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^ ' ABSTRACT 

The study of protoplanetary disks, where the planets are beheved to form, wiU certainly allow the formation of 
our Solar System to be understood. To conduct observations of these objects at the milli-arcsecond scale, infrared 
interferometry provides the right performances for T Tauri, FU Ori or Herbig Ae/Be stars. However, the only 
information obtained so far are scarce visibility measurements which are directly tested with models. With the 
outcome of recent interferometers, one can foresee obtaining images reconstructed independently of the models. In 
fact, several interferometers including IOTA and AMBER on the VLTI already provide the possibility to recombine 
' three telescopes at once and thus to obtain the data necessary to reconstruct images. 

In this paper, we describe the use of MIRA, an image reconstruction algorithm developed for optical interferometry 
data (squared visibilities and closure phases) by E. Thiebaut. We foresee also to use the spectral information given 
by AMBER data to constrain even better the reconstructed images. We describe the use of MIRA to reconstruct 
C3 , images of young stellar objects out of actual data, in particular the multiple system GW Orionis (IOTA, 2004), and 
discuss the encountered difficulties. 
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. 1. INTRODUCTION 

o ■ 

' The process of star formation triggered by the collapse and fragmentation of a molecular cloud is believed to lead 
— I to a protostar surrounded by an accreting disk undergoing important mass lost (wind, jets). The study of the 
' circumstellar disks during this pre-main sequence evolution is important because these disks are probably the birth 
OO , location of future planetary systems. Therefore, the understanding of their physical conditions is necessary before 
we can understand the planets formation process. The details of the physical processes which are taking place in 
^ , the circumstellar disk are not well known because of a lack of data to constrain them. The spatial distances at 
which the planetary, gas accretion and wind /jet formation occur, correspond to astronomical unit scales requiring 
milli-arcesecond spatial resolution to be probed unambiguously. Therefore, we used infrared interferometry which is 
' one of the techniques capable to reach this spatial resolution at those wavelengths. 

Interferometry consists of a coherent beam combination of the light received from an astronomical source and 
collected by two or more independent telescopes. As the light beams overlapped, we observe an interference pattern 
which is composed of fringes. An ideal interferometer measures the complex visibility V{u,v) which is, according to 
the Van Cittert-Zernicke theorem, the Fourier transform of the object brightness distribution I{x, y). In reality, the 
atmosphere turbulence induces random phase shifts and this phenomena does allow the Fourier phases to be properly 
calibrated. Therefore, we need to compute interferometric measurements that are insensitive to this perturbation : 
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• the squared visibility : V^^\u,v) = \ V{u,v) p 



• the bispectrum : B{ui, Uj,Vi, vj) = V(ui, Vi) x V(uj,Vj) x V(~Ui — uj, —Vi — vj). The phase of the bispeetrum 
is called the closure phase (CP), independant of the atmosphere delay phase. This method allows us to retrieve 
a part of the Fourier phase. 

From these interferometric data, the astronomers up to now are fitting simple parametric or homemade models in 
order to directly constrain their parameters to these visibilities and closure phases and in order to better understand 
the structure of the inner part of the disk. Yet, although important observations from these analyzes have lead to 
a profound revision of the inner disk structure at spatial scale corresponding to the planet formation, the lack of 
model independent imaging, without any knowledge of the observed structure, still forbids unambiguous diagnostics 
of the physical processes at work. 

With the outcome of interferometers able to combine 3 beams, like IONIC on IOTA and AMBER on the VLTI, 
and later of full imaging array like the second generation of instruments at VLTI, one foresees to obtain images 
reconstructed independently of the models. In fact, the AMBER instrument at the VLTI provides the possibility to 
recombine three telescopes at once and thus to obtain many data necessary to reconstruct images. 

2. METHOD 

In this section we present the method for performing image reconstruction. More details is given by Thiebaut^ in 
this volume. 

2.1 Image reconstruction in interferometry 

The observed bright distribution in a direction s' of an incoherent source is given by : 

/obs(s') = PSFix, y) * /t™c(a;, y) + N{x, y) (1) 

where PSF{x, y) is the point spread function, /true(a^, y) is the true object brightness distribution and N{x, y) is the 
noise; "*" denotes convolution. As interferometry does not measure in the image plane but in the Fourier space, we 
can compute the Fourier transform of the Eq. (1): 

Vohs{u, V) = 14rue(M, v) X S{u, v) + N' {u, v) (2) 

where S{u,v) is the Fourier transform of the point spread function, called the sampling function, Vtrue{u,v) is the 
Fourier transform of the true object brightness of distribution, called the true visibility and N'{u,v) is the noise in 
the Fourier space. 

The goal of the image reconstruction is allow to retrieve an image that can be trusted to be close to the real 
emission (/approx ~ /true)- To achieve this goal, the method is based on a Bayesian approach: to find a solution z 
which has the maximum probability given the data. Such solution is called the maximum "a posteriori" (MAP):^ 

z(^^P) =avgminf{z) (3) 



where the penalty function reads : 

f{z) = /data(2) + nfprior{z) 

This equation shows that we had to minimize a sum of two functions : 
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• the likelihood penalty term /data(z) which measures the compatibility of the solution z with the data. Typically, 
the data penalties are defined assuming a Gaussian statistics and it leads to a classical computation of a 
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• the prior penalty term /prior (2:) which allows us to account additional constraints not carried out by the data 
alone. In fact, there is more unknown parameters than measurements and an infinity of solutions are in 
concordance with the data. The regularization allows us to choice among all the compatible solutions the 
one which is the closest to an "a priori" . Because the noise usually contaminates the high frequencies, the 
smoothness prior is one of the most common regularization constraint: 



with ifc the discrete Fourier transform of x at fc-th discrete spatial frequency and where the weights w are 
chosen to ensure spectral smoothing to enforce smoothness of the Fourier Spectrum and thus compactness of 
the brightness distribution of the field of view. 

To control the importance of the priory penalty according to the likelihood penalty, a weight factor (see Eq. 4) 
is added in the equation. In practise, the best solution leads to fdata{z) — 1] where rj is the number of data. If 
fdata{z) is Small than ry, we should increase /x and if fdata{z) is larger than ry, we should decrease /x. 

2.2 Multi- Aperture Image Reconstruction algorithm (MIRA) 

The Multi- Aperture Image Reconstruction Algorithm (MIRA) is an image reconstruction algorithm devoted to opti- 
cal interferometry data developed by Thiebaut.^ Because of the small number of telescopes, optical interferometers 
are only able to measure very few Fourier frequencies at the same time and even less Fourier phases because of 
the atmosphere turbulence. The solution adopted in the MIRA algorithm is to directly fit the interferometric ob- 
servables (the squared visibilities and the closure phases) without explicity rebuilding the missing Fourier phases, 
unlike WISARD.^ The image reconstruction in MIRA is regularized by using various kind of regularization (entropy, 
quadratic and non-quadratic smoothness ...). MIRA states the image reconstruction as an inverse problem solved 
by minimizing a so-called penalty function (see Sect. 2.1). Owing to the specific relationship between the object 
brightness distribution and the interferometric data, the penalty to minimize cannot be guaranteed to be convex. 
At this time, the algorithm for the optimisation used in MIRA is the VMLM-B non-linear constrained optimization 
algorithm^ and which imposes positivity and normalization of the sought image but only yields a local minimum. 
Hence the initial image of the iterative algorithm is a deciding parameter of the algorithm. 
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3. IMAGE RECONSTRUCTION ON YOUNG STELLAR OBJECTS 



In this section, we present the results obtained on two objects. The first one is a known young binary, GW Orionis, 
and the other one is a young star which is believed to have disk, HD 45677. 



Table 1. Characteristics of the classical T Tauri star GW Orionis 



RA (J2000) Dec (J2000) my m./ mg xhk Dist. (pc) Spec, type 
5h39m08s +11 °52'12.7" 9J TJ 7l 6^6 300 K3V 



3.1 GW Orionis 

3.1.1 Characteristics of GW Ori 

GW Orionis (see characteristics on Table 1) is a young stellar object and belongs to the classical T Tauri group. 
Mathieu et al.^ discovered that this star was in fact a spectroscopic binary with an orbital period of 242 days. 
The masses of the binary component is 2.5 solarmasses for the primary and between 0.5 and 1 solarmasses for the 
secondary and the separation is slightly more than 1 AU. The ccnter-of-mass velocity has been observed to vary 
over a period of 1000 days, suggesting the presence of a third star in the system or an m = 1 perturbation in the 
associated disk. 

The observed spe(;tral energy distribution (SED) of GW Orionis shows a large near- and far-infrared excess over 
the stellar photosphere and presents two characteristics : 

1. the morphology of the SED is double peaked with a minimum in the continuum near 10 /xm; 

2. there is a very strong 10 fivci silicate emission feature. 

The interferometric data on GW Orionis used for the image reconstruction (see Fig. 2 in red squares for the 
data and Fig. 1 on the left panel for the [u, v) plane) have been taken at the IOTA interferometer.^ Located on 
Mount Hopkins in Arizona, IOTA is a long baseline interferometer that observes at visible and NIR wavelengths. It 
has three telescopes movable among 17 stations along two orthogonal linear arms. By observing a target in many 
different array configurations, IOTA can synthesize an aperture of 35 x 15 m (corresponding to an angular resolution 
of 5 X 12 mas at 1.65 /^m). 

3.1.2 Reconstructed image with MIRA 

The image reconstruction of GW Ori has been done with the following parameters: 

• a pixel size of 0.1 mas to have a more precise image 

• an image width of 400 pixels to have a sufficient field of view 

• a regularization of type "smoothness" 

The result of the image reconstruction is shown in Fig. 1. Besides the image, MIRA gives also a comparison 
between the measured interferometric data (squared visibilities and closure phases) and the reconstructed interfer- 
ometric data computed by the Fourier transformation of the image, which can give us an approximation of the 
and show which data point are not matched by the image reconstruction (see Fig. 2 in green). 

We find that the GW Ori system appears to be a triple system, with a binary already known (see section above) 
and a faint third component. Prom this image, we can measure the different parameters of the triple system, namely 
the flux ratio between a reference (taken as the brightest one) and the two other components, the distance between 
the components and the position angle. The parameters measured on the image are shown in line 1 of Table 2. 
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Figure 1. («, v) coverage of the data set and the result of the reconstructed image by MIRA for GW Ori (with the dirty map). 



Table 2. Results of the different methods used on the GW Ori data, 1 indicates the reference component, 1-2 the thight binary 
and 1-3 the wide binary. 
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Figure 2. Graphical result of the different methods used on GW Orionis. Left : the squared visibilities; Right : the closure 
phases. In red squares, the interferometric data; in green triangles, the reconstructed interferometric observables from the 
MIRA reconstructed image; in blue circles, the binary fit; in black stars, the triple point fit. 



3.1.3 Independant model fitting 

In order to evaluate our confidence on the image reconstruction, we decided to perform a model fitting totally 
independant of MIRA results and to convince ourselves that the both methods lead to the same results. We followed 
a two-step fitting procedure : 



1. Binary fit : in a first step we carry on a grid search of the minimum x^- For that purpose we compute an 
hypercube of minimum exploring all realistic binary parameters. From the resulting hypersurface, we 
determine an approximate position for the global minimum (see Table 2, line 2 and Fig. 2 in blue circles). The 
graphical representations of the cube (see Fig. 3) allow us to evaluate the sharpness of the peak and 
the possible occurence of local minima; the curvature of the chi-2 around the minimum also provide estimation 
of the precision of the parameters. As shown Table 2 (line 2), the of this method is not satisfying. 

2. Triple point fit : in order to improve the x^ and motivated by the spectroscopic clear companion evidence, we 
add a third component to the model. This second step consists in a classical Levenberg-Marquardt minimisation 
of the x^ with the six free parameters, three for each binary considering one component as the reference. The 
result of the binary fit is used as initial guess for the wide binary's parameters, whereas the parameters of the 
tight binary are choosen randomly. The LM fit result is shown Table 2 (line 3) and graphically Fig. 2 (in black 
stars) . 



We see that the result for the wide binary is approximately the same with both methods, whereas the results for 
the tight binary are quite different. The difference can be explained by the limited resolution of the interferometer 
and also by the fact that our data constrain most the long separation than the shorter one due to the longer of 
baselines used for the observations. 
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Figure 3. Graphical representations of the cube between the GW Ori data and a binary model for a whole range of values 
for the three parameters : the flux ratio, distance and position angle. 

Table 3. Characteristics of the Be star HD 45677. 



RA (J2000) Dec (J2000) my mj mg Dist. (pc) Spec, type 

6h28ml7.4s -13°3'11.1" %m 7^2 O 48 350 Bpshe 



3.2 HD 45677 

HD 45677 is a YSO that belongs to the Ae/Be Herbig group (see its characteristics Table 3). The first data of 
HD 45677 have been published by Monnier et al/ given the first closure phase survey of YSOs using an infrared 
interferometer. In this paper, they did parametric imaging on HD 45677, proposing this Be star to be surrounded 
by an elongated and highly skewed dust ring (see then right panel of Fig. 4) . 

The interferometric data of HD 45677 used for the reconstruction are those of Monnier et al.^ (see Fig. 5 in 
red for the data and Fig. 4 for the (u, v) plane) which have been obtained also on the IOTA interferometer, as GW 
Orionis (see Sect. 3.1) 

The image reconstruction of HD 45677 has been done with the parameters below : 

• a pixel size of 0.2 mas 

• an image width of 256 pixels 
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Figure 4. (u, v) coverage of the data set and the parametric imaging of IOTA data for HD 456777 
• a regularization of type "smoothness" 

The result are showed Fig. 5. The image seems to show some points of a skewed disk, agreeing with the results 
obtain by Monnier et alJ To obtain a more precise image, we need more interferometric measurements, especially at 
longer baselines. This is why we requested several nights of observations at VLTI with different triplets of baselines 
in order to have a good coverage of the (m, v) plan. 

4. CONCLUSION AND FUTURE 

With only a few algorithms in the world, the image reconstruction in optical interferometry is really at its first steps 
especially for YSOs. The problem is much more complex than the image reconstruction in radio interferometry: in 
fact, in optical interferometry, we do not measure the complex visibility directly, the presence of the atmosphere 
leads to a lack of information on the phase and finally we have less measurements at the same time. Therefore, 
the optical interferometry needs its own method to retrieve images from the interferometric observables (squared 
visibilities and closure phases). 

With one of these algorithms, MIRA, we show in this paper the first images from real long-baselines interfero- 
metric data of the environment of young stellar objects. The results have been analyzed with a new approach: we 
foung good correspondence between the reconstructed image and the parametric images, especially for GW Orio- 
nis. This demonstrates that this technique allows us to raise the degenerancy which appears when we analyse the 
interferometric data with model fitting. However, as the image reconstruction is too heavy to do statistics on the 
reconstructed image and as there is not a method yet which can compute a significant error in order to estimate the 
rightness of the reconstructed image, we have difficulties to easily derive scientifical results from the reconstructed 
image. 

The image reconstruction in optical interferometry needs therefore to make progress but it already shows that the 
technique will give great facilities in the analysis of the interferometric data. With future instrument as the VLTI 
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Figure 5. The result of the reconstructed image by MIRA for HD 45677 (with the dirty map) . Second hne left : in red squares, 
the measured squared visibilities; in green triangles, the squared visibilities of the reconstructed image. Second line right : 
idem for the closure phases. 



Spoctro Imager VSI, a second generation instrument for the VLTI, we can make image reconstruction of several 
astronomical objects® : with the possibility of observing with 4 to 6 telescopes and with a groovy decrease of the 
error bars, the {u, v) plane coverage is well much better and the reconstruction made by MIRA is easier (see Fig. 6 
for an example on YSOs) 
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Figure 6. Reconstructed image from simulated data for the science case of the VLTI Spectro-Imager (VSI): YSOs. First hne: 
the {u, v) coverage of the data set. Second hne: the simulated image (left, logarithmic scale) and the reconstructed image 
by MIRA (right). Third line: the squared visibilities (left) and the closure phases (right) in red squared for the data and in 
green triangle for the MIRA results. 



